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Abstract 

We study the change of phase space structure of the rectilinear three-body problem when the 
mass combination is changed. Generally, periodic orbits bifurcate from the stable Schubart periodic 
orbit and move radially outward. Among these periodic orbits there are dominant periodic orbits 
having rotation number (n — 2)/n with n > 3. We find that the number of dominant periodic orbits 
is two when n is odd and four when n is even. Dominant periodic orbits have large stable regions 
in and out of the stability region of the Schubart orbit (Schubart region), and so they determine the 
size of the Schubart region and influence the structure of the Poincare section out of the Schubart 
region. Indeed, with the movement of the dominant periodic orbits, part of complicated structure of 
the Poincare section follow these orbits. We find stable periodic orbits which do not bifurcate from 
the Schubart orbit. 

1 Introduction 

In the present study, we continue our work on the rectilinear three-body problem (Saito & Tanikawa, 
2007; hereafter referred to as Paper I). The rectilinear three-body system is such that three particles are 
on a line. This problem has two degrees of freedom. 

The structure on the surface of section of the rectilinear three-body system has been studied using the 
Poincare map (Hietarinta & Mikkola, 1993; hereafter HM1993). They found that the Poincare section is 
divided into three basic regions: the Schubart region, the chaotic scattering region, and the immediate 
escape region. Inside the scattering region, the interplay time is so sensitive to the initial conditions that 
there seemed to be no structure. HM1993 also found that the number of scallops, which constitute the 
immediate escape region, increases as the mass of the central particle becomes smaller. 

Tanikawa and Mikkola (2000a; hereafter referred to as TM2000a) studied the structure of the Poincare 
section for the equal-mass case using symbol sequences which record the collisional history of orbits. They 
used the Poincare section as the initial condition surface and associated with each point of the surface 
its future history including the final motion. They then succeeded in showing that the chaotic scattering 
region is filled with points whose orbits end with triple collision and that these points form well stratified 
curves. TM2000a's fine structure of the Poincare section for the equal-mass case shows that (1) stratified 
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triple collision curves make a sector together with the sub- region of the immediate escape region (i.e., 
scallop); (2) four sectors surround the Schubart region. 

The study of a similar dynamical system, the collinear Coulomb three-body problem with electron- 
ion-electron configuration, carried out by Sano (2003) is helpful to understand our system. In his system, 
there exist the Schubart region and unstable periodic points on its vertices. He confirmed the existence of 
the separatrices of these unstable points in his system numerically following the mapped points starting 
near the vertices. 

In Paper I, we studied how the structure on the Poincare section changes as the central mass varies. 
The chaotic scattering region is foliated with triple collision curves in general mass combination. This 
foliation makes two types of block: arch- and germ-shaped blocks (hereafter simply arches and germs; 
see Fig. []}. The chaotic scattering region consists only of arches near the mass combinations of total 
degeneracy, whereas it consists of arches and germs for the mass combinations away from total degeneracy. 
Here by total degeneracy we mean that stable and unstable manifolds of the two fixed points on the Triple 
Collision Manifold connect smoothly (McGehee, 1974). In the totally degenerate cases, the solutions can 
be analytically continued beyond triple collision. The number of arches always coincides with the number 
of the scallops. The mass combinations of total degeneracy are numerically found by Simo (1980). We 
introduced a family of sets of symbol sequences, and partitioned the Poincare section by sets of points 
corresponding to sets of symbol sequences. We then found that an arch consists of sets of points whose 
symbol sequences are arranged according to a simple rule (see §4.1 of Paper I). The composition of 
arches changes each time the number of scallops increases, that is, the mass parameters pass through 
a totally degenerate case. We observed that a number of germs bifurcate from arches and then germs 
from different arches construct new arches. However, we could not understand the dynamics behind the 
behaviour of germs. 

In the present paper, we also study the structure change of the Poincare section as the central mass 
varies (we only consider the symmetric mass combinations). However, in contrast to Paper I, where 
the structure of the Poincare section was considered in relation to triple collision, here the structure is 
analysed in relation to the motion of periodic orbits bifurcated from the Schubart orbit. From another 
point of view, the azimuthal structure of the Poincare section with centre at the Schubart point has been 
studied in Paper I, whereas the radial structure with centre at the Schubart point will be studied in the 
present paper. 

The rotation number with respect to the (Schubart) fixed point will be introduced. General tendency 
of the motion of the bifurcated orbits from the fixed point is analysed. The rotation number monotonically 
decreases as a increases. As we shall see, there is the sequence of rotation numbers such that the 
corresponding periodic orbits has significant influence on the Poincare section, these orbits will be called 
dominant periodic orbits. We mainly concern with these orbits. As the mass parameters change, the 
dominant periodic orbits bifurcate, recede from the fixed point, and their stability region shrink. We 
will relate these behaviours of periodic points to the growth of the germs observed in Paper I. 

The organisation of the present paper is follows. In §2, we introduce the equations of motion of our 
system. We define the Poincare section and introduce the rotation number in the fixed-point-centric 
coordinates. The §3 is for the results. In §3.1, the bifurcation of periodic points is analysed as the mass 
parameter changes. In §3.2, we find that periodic points with special rotation number are dominant over 
the structure of the Poincare section. In §3.3, we show in detail the radial motion of the periodic point 
with the growth of the germs as the mass parameter changes. Finally in §4, we summarise our results. 
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2 Method 



2.1 Equations of Motion and Initial Points 

We introduce the equations of motion, the Poincare section, and parameters describing mass combination 
of the particles (MH1989). Then, we introduce symbol sequences recording the collisional history of orbits 
(TM2000a). 

Let mi, mo, and TO2 denote the masses of the three particles from left to right on the line. Let 
(gi, g2) be the mutual distances to-ito-o and TO0TO2, and {pi,p 2 ) the conjugate momenta to (gi, g2). The 
Hamiltonian H is written, with the kinetic energy K and the force function U , as 

H = K-U, (1) 
2 Vtoi mo/ 2 Vtoo TO2 / too 

TOiTOq TOqTO.2 TO1TO2 



u 



qi q_2 qi + qi 



In this paper, we consider the case that total energy E, where H = E = const., is negative. We restrict 
ourselves to E = — 1, since negative energy system can be brought to E = — 1 using the homogeneity of 
U. We here adopt variables used in MH1989. For readers' convenience, we repeat the formulation. Let 
us consider the canonical transformation from (qi,pi,t) to (Qi,Pi,f) defined by 

qt = Ql p i = 2p i Q i (i = l,2) (2) 
dt — qiq-2dt' . (3) 



Then the new Hamiltonian reads 



T = qi q 2 (H-E), (4) 



and the equations of motion become 

dt> ~ dP^ df ~ dQ t w 

The transformed equations (O show that binary collisions (Q\ = 0, Q2 7^ or Q2 — 0, Q\ 7^ 0) are 
regularised, whereas triple collision (Qi = Q2 — 0) is still singular. The solution to Eq. ([5]) is continued 
beyond binary collision. A binary collision is interpreted as an elastic bounce in the physical coordinates. 

The equations of motion allow a solution qi(t) = Aga(t) for all t for which the solution is defined. 
This is called the homothetic solution (Irigo yen and N ahon, 1972). Here A is a constant depending on 
masses. This dependence is given by the following equation: 

—mi + ztoo + TO2 z 5 — 2z 3 + I7z 1 + z , . 

where A = . (6) 



TOl + TOo + TO-2 z 4 — 10z 2 — 7 1 

The intersection of the energy hyper-surface H(q,p) = E and the hyper-surface gi = Ag2 is a two- 
dimensional surface. We take this surface as the Poincare section n. In the following, we introduce 
coordinates (9, R) on n. The variable R is defined by 

R = -(gi + g 2 ) where gi = Ag 2 . (7) 

With solving Eq. ([7]) to obtain gi and g2 for a given R, the value of K is determined as 

2Ai? 2R 

K = E + U{ qi ,q 2 ) with gi = y - ^' 92 = YTA' (8) 
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For this value of K, the parameter 9 determines the ratio of q\ to q 2 : 



VK cos 9 = \J A — Kr 2 2 qi — \J B — nri 2 q 2 

\Z~Ksm9 = r 2 \ r nqi - r\sJ~Kq?,, (9) 



where 



m 1 (roo+ro 2 ) D m, 2 (m + mi) m 2 m 

A= 2M ' B= 2M ' C= W M = m 1 +m +m„ 

2A 2 4AB-C* 2 

r l = 1 7; r 2 = 1 7, « = ~TT~a o ™ o ^ 7; < 6> < 27T. 

A + l A + l 4(Ari 2 + Br 2 2 + Cr Y r 2 ) ~ 

The variable R takes its maximum value i? max when K = Q. Substituting if = in Eq. l[8"]). we 
obtain 



R -M 



A + l , 

toottii + (A + l)mom2 + m\m 2 . (10) 



A 



We denote by II the side of the Poincare section with < 9 < tt and by II* the other side with ir < 9 < 2ir. 
By the definition of 9 in Eq. q(9 + tt) = —q{ff), where we denote by q(9) the (91,32) for given 9. 
Since the original equations of motion are invariant under the transformation (<&,£) — » (— q%,— t), it is 
enough to integrate orbits toward the future for each point (9, R) in II U II*. Nonetheless, we only study 
the structure of the surface II for future integration. The reason why we do so is that an orbit intersects 
alternately with II and II* until the corresponding triple system disintegrates into a binary and a single 
particle. Only certain orbits starting from the immediate escape region intersect with only II or II*. 
Since we are interested in the structure of chaotic scattering region, such escape orbits can be neglected. 

We assume the total mass of the particles to be three without loss of generality and introduce 
parameters a and b to represent the masses: 

mi = 1 — a — 6, mo = 1 + 2a, m 2 = 1 — a + 6, (11) 
where a > -1/2, b > 0, a + b < 1. 

The parameters move in a triangular area of the (a, 6)-plane. We will call this a mass triangle. In Paper 
I, we saw the structure change of the Poincare section over the mass triangle. In the present paper, we 
only consider the symmetric mass configurations: 6 = 0. The reason is that we saw in Paper I that 
asymmetry of the mass configuration in many cases does not add topologically new structure to the 
symmetric case. 

There are three types of collision in our system: left-centre (gi = 0, q 2 7^ 0), centre-right (q 2 = 
0, q\ 7^ 0), and triple (gi = q 2 — 0) collisions. Let us denote these collisions by symbols '1', '2', and 
'0', respectively. An orbit of the partly regularised equations Eq. ((SJ) repeats binary collisions until a 
triple collision takes place. These collisions can be recorded as a sequence of the symbols. We study 
the evolution of orbits using the symbol sequences instead of orbits themselves. It is to be noted that 
we encode into symbol sequences the future behaviour of the orbits starting at points of the Poincare 
section. So symbol sequences are not bi-infinite but singly infinite to the future. 

Integration of orbits are carried out for grid points. We introduce the grid points (9i, Ri) with different 
resolution on II: 

(9i, Ri) = (tt x - — ,R max x — — ), where i = 0, • • • , 9 div - 1, j = !,-■■, R div , 



^div Rdr 

and / (i) 6> div = 540, i? d iv = 300 

(ii) 9 div = 1800, i?div = 1000 



(12) 
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X. Foot- points 



Figure 1: The schematic illustration of the structure on the Poincare section. The Poincare section 
consists of the Schubart region (I) and sectors (V). A sector (V) includes one scallop (IV), one arch (VI), 
and optional germs (VII) and pipelines (VIII). In the third sector from the left, it is shown that both 
an arch and a germ consist of the strata of triple collision curves. In the blank area of the third sector, 
there are other germs and/or pipelines, although they are not drawn. The pipeline (VIII) is an object 
which connects two germs in different sectors. Note that there are narrow escape regions also inside the 
sectors. 



Most of structures on the Poincare section and their change due to the mass variation can be resolved with 
resolution (i). Detailed structures may be lost, if one uses resolution (i). However, once we understand 
certain fine structures with resolution (ii), we can use resolution (i) to follow global behaviour of the 
structure of the Poincare section when mass parameters are changed. We have obtained the first 64 digits 
of symbol sequences for the grid points through integration of orbits. These constitute the basic data in 
the rest of the paper for the distribution of symbol sequences on II. We integrate Eq. ([5]) using DIFSY1, 
which is an implementation ( Press and T cukolsk y, 1999| of Bulirsch-Stoer method. In our preliminary 



research ( |Saito fc T anikawa, 2004), we introduced cylinders (sets of symbol sequences with given words 



where a word is a finite sequence of symbols) S c j, and S, 

f {(21)«(2)J -...|i> 0, j > 1} forc = 2 J + l 
° c ' j ~ \ {(21) i (l^'---|i > l,j > 1} forc = 2i 1 ^ 

= / {(21) i (2)~|i>0} forc = 2i + l 
c '°° ~ \ {(21) l (l)°°|i > 1} forc = 2i [ ' 

S c — ^-ij<ooS c ,j ■ (15) 

Let reg(c) and reg((c, oo)) denote the regions on IT whose symbol sequences belong to S c and S CiOC . The 
three basic regions in HM1993 correspond to our regions as follows: 

Immediate escape region reg((c, oo)) 1 < c < npop 

Chaotic scattering region reg(c) U reg((c', oo)) 1 < c < oo, npop + 1 < c' < oo . 

Schubart region reg(c) c = oo 

For later convenience, we recall the structure of the Poincare section and names for its elements 
obtained in Paper I (Fig. [T]). The separatrices running from the vertices of (I) the Schubart region divide 
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the outside of the Schubart region into (V) sectors. A sector includes (IV) one scallop. In the remainder 
of the sector, (IX) triple collision curves are foliated. This foliation makes one or more blocks. One block 
is arch-shaped (VI, arch) and the other blocks are germ-shaped (VII, germ). A germ bifurcates from an 
arch, as a changes. Germs that bifurcate from different arches finally gather to construct a new arch and 
to reconstruct other arches. As is shown later, the growing process of germs is related to the movement 
of the periodic orbits bifurcated from the Schubart orbit. 

2.2 Poincare map and rotation number 

We here define the Poincare map on the Poincare section (9,R). An orbit starting from the Poincare 
section repeats the intersection with II and II* alternately. When an orbit intersects with II (resp. II*) 
at {6,R) and II (resp. IP) again at (0",R"), we define a map T from (9,R) to (0",R"). 

T : {9, R) ^ {9", R"), where (9, R), (9", R") G II or (9, R), (9", R") e II* (16) 

Equation {5} has a simple periodic solution, the so-called Schubart orbit (Schubart, 1956), whose 
symbol sequence is (21)°°. The intersection Pq = {6q, Rq) of the Schubart orbit with II is a fixed point of 
the map, namely Pq = T(Pq). The linear stability of T around Pq over the mass triangle is studied by 
HM1993. Approximately in the mass combinations such that m\ < tuq < m2, Pq is hyperbolic whereas 
in the other combinations it is elliptic. The stability region is called the Schubart region. A point P 
such that P — T q (P) is called a q-periodic point and the sequence {P,T 1 (P), ■ ■ ■ ,T q ^ 1 (P)} is called a 
q-periodic orbit. We are interested in the periodic orbits that bifurcate from the fixed point as the mass 
parameter a is changed. 

In order to describe the elliptic motion around Pq, we introduce the rotation number. The rotation 
number is the averaged number of rotations per iterate of T. First, we introduce polar coordinates (£>, A) 
with centre at Pq = (9q,R ): 

DcosA = (9-9o)g e , DsmA = (R-R Q )g R , 



where gg — 180/-7T and g R = 100/i? max . The values of gg and ga are arbitrarily chosen such that 
both 9 and R equally contribute to the values of D and A. Now the Poincare map is expressed as 
T : (D,A) i— > (D',A'). We then introduce the effective rotation number v n (D,A;a) and the (exact) 
rotation number ^^(D, A; a) at (D, A] a) in the following equations: 

diffU' A)-i A '~ A (if A' - A > 0) 

1 n 

v n {(D 1 A)-a) = —Y j m{A^,A^), where A«) = T\D, A) (18) 
i=i 

Vao((D i A)' ) a)= lim i/k(D,A;a) (if the limit exists). (19) 

k — >oo 

We here remark several points related to the numerical calculation of the rotation number. The first 
point is the selection of n of the z/„((Z?, A); a). In our approach, we first take an fixed integer N and 
iterate the mapping N times. We then find a number n (1 < n < N) such that diff ( A^ , A) = 
mini<i<Ar{diff {A®, A)}. The second point concerns the rotation number of periodic points. Suppose 
p and q are integers. If P is a periodic point, then there exist integers p and q such that v q (P;a) — 
^oo(P; a) = p/q (by the above method of selection of n, n — q). This coincidence of the effective and exact 
rotation numbers helps us to find the periodic points: the points on the contour with v q (P\ a) = p/q on 
the the Poincare section are the candidates of the periodic points with ^(P; a). The third point is that 
the rotation number v oq {Pq]o) is calculated from the eigenvalues of linearised map of T aX Pq. When 
this value is rational, the fixed point and the periodic points are degenerated. If the mass parameter a 
is changed, the periodic points bifurcate from fixed point Pq- 
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Figure 2: A Poincare map in order to select candidates for periodic points. In this map, a = —0.25 and 
the expected periodic points are those with v — 4/6. The centres of elliptic motions under the Poincare 
map are the candidates. 

2.3 Search for the periodic points 

We explain how to detect the periodic points. First, we have to calculate lim/j^o Voo{{D, A); a) for 
various values of a. For a given rotation number p/q, we can determine the mass parameter a where 
corresponding periodic points appear: solve Voo{Pq] a) = p/q for a. Second, we find candidates at a mass 
parameter a little distant from the exact bifurcation. The candidates for periodic points are found from 
the Poincare map such as is shown in Fig. [U The candidates for elliptic points can be taken at the 
centre of the libration and for hyperbolic points at the saddle between two librations. Finally, we obtain 
a periodic point through the Newton-Raphson method Eq. (|20|) starting from an arbitrary point on the 
contour. This method is used in MH1991 to find the fixed point. 

X n +\ = X n — 

3 Results 

We now show the behaviour and influence of the periodic points on the Poincare section. First, their 
bifurcation and radial movement away from the Schubart orbit will be discussed. The values of parameter 
a where the periodic points bifurcate can be determined from the rotation number, ^^(Pq; a), of the 
Schubart orbit as a function of a. The relation between radial movement and the variation of a will 
be studied. Second, the number of orbits which appear in respective bifurcations for special rotation 
numbers will be discussed. Finally, the influence of such special periodic orbits on the Poincare section 
will be discussed. 

3.1 Bifurcation as the parameter changes 

In this subsection, we show how periodic points bifurcate and move outward from the fixed point Pq, 
when the parameter a is changed. As is stated in section 2.2, period-g points bifurcate if Voc{Po', a) = p/q 
with integers p and q. In order to know when the periodic orbits bifurcate, we calculate v co {Pq]o) as a 
function of a (Fig. [3]). For several rotation numbers, we show the values of a in Table[TJ Figure [3] shows 




dT(X n ) 

dX n 



- 1 



(T p (X n ) - X n ) = 0. 



(20) 
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Vao(P{)]a) is monotonic: 



V 



dvoo (Pa; a) ^ Q 
da 

rational number p/q, the periodic points with = p/q bifurcate 
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(21) 
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Figure 3: The rotation number v(Pq 1 a) as a function of a. 



all. Note that this graph shows that the period-3 orbit bifurcates at a = ax/3 —■ 0.41. According to 



^00 1/4 


1/3 


2/4 


3/5 


4/6 


3/7 


a 0.61 


0.41 


-0.019 


-0.15 


-0.25 


-0.31 



Table 1: The values of a at which the periodic points bifurcate for several rotation numbers 



Moser (1958), generally in Hamiltonian systems when the period-3 orbit bifurcates, the mother periodic 
orbit becomes unstable. As an application of this result, he tried to explain the origin of the Kirkwood 
gaps. Also in our system, the instability at the bifurcation of the period-3 orbit was confirmed by 
Hietarinta and Mikkola (1993). We have followed the position of the periodic points with ia^ = 1/3 on 
the Poincare section with continuously changing a. As a result, when a is decreased, with a passage 
a = ax/3, the triangle formed by these periodic points approaches (a > ax/3), collides (a = ax/3) with, 
and recedes (a < ax/3) from P$ (Fig. Ufa)). The orientation of the triangle is inverted before and 
after ax/3 (Fig- 1Kb)). Similar behaviour was observed by Henon (1970) in Hill's case of the restricted 
three-body problem. In his case, mother periodic orbit is the retrograde satellite orbit. 

We want to know the direction of change of a for the bifurcated periodic points to move outward 
from Pq . Suppose that the periodic points bifurcate at a = a p / q and exist at a = a p / q + da as the points 
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(a) The distance from Pq 
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(b) The inversion of triangle's orientation 




Figure 4: The radial movement of the periodic points with = 1/3 from Pq. (a) the distance between 
one of these point and Pq. 



with a finite distance dD from Pq. The sign of da can be obtained from Eqs. pTj) and (|23|) . as follows. 
The total differential of v as a function of a and D at Pq is 



dv = 



(-) 

\da) 



P da+ (^) 



dD. 



Since dv = on periodic points, we obtain 

da = 



dv / dv 
~dDl ~da~ 



dD. 



(22) 



We only need the sign of dv 00 (PQ-,a)/dD, since we already know the sign of dv O0 (PQ-,a)/da. We 
approximate the derivative dv n (P;a)/dD by the difference (v n ((D + AD, A); a) — v n ((D, A); a))/ AD 
for small AD. Since dv n (PQ-,a)/dD should be independent to A, we select arbitrarily the value of A: 
A = 0. Moreover, we have Voo(Pq; a) '= v n (PQ-, a) for large n < N: N — 256. Finally we put AD = 0.1. 
Figure [5] shows 8v 00 (Pq; a)/dD as a function of mass parameter a. The figure shows that the sign of 
dv oc (P ;a)/dD is 



duoo(Po;a) 
dD 

dD 



> (for a > aif 3 ) 
< (for a < ai/ 3 ) 



(23) 



Taking into account Eqs. ipTjl and (|23p . and dD > 0, we find that da is positive for a > a 1 / 3 and 
negative for a < ai/ 3 . The recession from a — a l / 3 , irrespective of the sign of a, yields successive 
recessions of the periodic points from the Schubart orbit. This implies that absense of the periodic 
points at a — ai/ 3 (except for the periodic points independent of Pq). 



3.2 Dominant Periodic Orbits 

In the present subsection, we will find periodic orbits and follow their motions for various rotation 
number. These results show us that the periodic orbits with rotation number [n — 2)/n dominate the 
structure of the Poincare section. We call these the dominant periodic orbits. We discuss a few features 
of the dominant periodic orbits. On the other hand, periodic orbits with the other rotation numbers have 
too small stable regions to numerically find the location on the Poincare section, so we do not consider 
these. 
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Figure 5: Dependence of dv(D : a)/dD on a. The sign of dv(D, a)/dD changes at a 



H/3- 



According to our numerical results, there is a rule for the number of dominant periodic orbits. Suppose 
that the rotation number of these orbits is = (n — 2)/n. Indeed, for even n, the period is not n but 
n/2, and therefore, if we would strictly write the rotation number, it should be = {n/2 — l)/(n/2). 
We have found the following rule for the number of periodic points. 



A Rule For the number of dominant periodic orbits with = (n—2)/n, the following rule is confirmed 
from n — 4 to n = 19. 

• the number of orbits is two for odd n. 

• the number of orbits is four for even n. 

We show an example of this rule in Fig. E] (n = 5 and 6) and Fig. \7\ (n = 17 and 18). For the case of 
fi = 3 (iSco = 1/3), only one unstable orbit appears through bifurcation as already seen in Fig. [UJb). It 
is important to note that according to the Poincarc-Birkhoff theorem (Birkhoff, 1913) the number of the 
periodic orbits is 21 with an integer £, and in the case of the standard map ( Chirikov, 1979 ) £ is believed 
to be one. We do not know the number of non-dominant orbits. 



(a) v = 3/5 



a =a/B PPs (a=-0.1B) 





Figure 6: Periodic orbits for is^ = (n — 2)/n with (a) n = 5 and (b) n — 6. These show the numbers of 
orbits are 2 for odd n and 4 for even n. We connect the periodic points which belong to the same orbit. 
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(c) v= 15/17, (a, b) = (-0.461844,0) 



(d) v= 16/18, (a, b) = (-0.466010,0) 




Figure 7: Periodic orbits for = (n—2)/n with (a) n = 17 and n = 18. For each orbit, its orbital points 
are plotted by different symbols: single and double circles (stable orbits), multiplicative and additive 
symbols (unstable orbits). Two symbols appear for odd case (n = 17), whereas four symbols do for even 
case (n = 18). The background dots shows the Poincare map of segments of several lines running from 
Po- 

3.3 Periodic orbits and the structure on the Poincare section 

In the present subsection, we study how the dominant periodic orbits are related to the structures on 
the Poincare section. 

Because of the reversal of the bifurcation at a = <2]/3, we observe different movements of periodic 
orbits for two cases: da > for a > ax/3, an( i da < for a < ai/ 3 . In the latter case, there is a nearly 
periodic change of the Poincare section, and one cycle is from the bifurcation with = (n — 2)/n to 
that with = {(n + 1) — 2}/(ra + 1). The change in the former case is basically similar to the one 
cycle in the latter case (Note that for a > ai/3 only = 1/3 has the form = (n — 2)/n). Hence, we 
propose the scenario for the latter case only. The following scenario is based on the observation of the 
cycles starting respectively at Voo = 2/4,3/5, 4/6, 5/7, and 6/8. 

We explain our scenario referring to Fig. [S] which corresponds to the cycle starting at = 3/5. 
The bifurcated periodic points with is^ = (n — 2)/n (n > 4) recede from Po as a decreases. As the 
distances become large, the stable regions around the periodic points, called islands, influence the shape 
of the Schubart region. In both Fig. [5] (a) and (b), five triangular islands can be seen. The separatrices 
which connect the unstable periodic points (and envelope the stability region of the stable ones) is of 
polygramic shape. The germs (see Section 2.1) grow along the separatrix. That makes the Schubart 
region be of polygramic shape. When the periodic points together with their islands get out of the 
Schubart region, the Schubart region returns to the polygonal shape. The germs now intrude between 
the Schubart region and islands. As a further decreases, the stable periodic points sink into the gaps 
between arches, and their islands shrink. The germs grow and follow the sinking stable periodic points. 



11 



As a result, a set of germs pile up over an arch. The piled germs become a new arch. On the contrary, the 
unstable periodic points stay around the vertices of the Schubart region: their separatrix approximates 
the boundary of the Schubart region. If we continue to decrease a, the next dominant periodic orbits with 
^oo = {(n + 1) —2}/(n + l) bifurcate. The evolution of the structure of the Poincare section following this 
bifurcation is similar to the case v = (n — 2)/n. It is to be note that there seems to be a correspondence 
between the structure around the Schubart orbit (the fixed point) and a structure around a bifurcated 
periodic point. Around the Schubart orbit, there are the Schubart region and a number of arches, while 
around a bifurcated periodic point there are the islands and the region filled with germs. 

In what follows, let us see numerical evidence that supports our scenario. For this purpose, we follow 
the change of the Poincare section over the range of a corresponding to ia^ =3/5 and 4/6. In Fig. 
the periodic points on the Poincare section for the range —0.150 > a > —0.250 are displayed. The 
structure of the Poincare section is represented in two ways: the triple collision curve and the partitioning 
according to the cylinders defined by Eq. (fl~5"|) . The separatrices of the periodic points are also plotted. 

The radial recession of the periodic points and the enlargement of the islands, as a decreases, can 
be seen in Fig. [8] (a) a = —0.150 and (b) a = —0.160. The germs grow along the separatrix which 
encloses the islands at (b) a = —0.160, whereas the germs split the Schubart region and the islands at 
(c) a — —0.166. This means that the periodic points get out of the Schubart region at a' such that 
—0.166 < a' < —0.160. After getting out of the Schubart region, islands become smaller. This is seen 
from the transition from (c) a = —0.166 to (d) a — —0.170. 

The comparison of (d) a = —0.170 and (e) a = —0.186 shows us that the germs grow and follow 
the stable periodic orbit which is sinking between arches, and this growth of the germs results in the 
reconstruction of arches. For example, the second arch from the left has reg(2), reg(6), and other small 
regions at (b) a = —0.160, whereas reg(2), reg(7), reg(12), and other regions at (e) a = —0.186. Note 
that the head of the germ reg(5, • • •) (the upper right in (e)) does not reach the 0-axis and the number 
of arches is still four. The reconstruction process was studied well in Paper I. However, we have now 
found that this is promoted by the periodic orbits. Again in (e) a — —0.186, the unstable periodic points 
marked as five open circles stay at the vertices of the Schubart region. On the other hand, the stable 
periodic points marked as remaining five open circles are on the midway to the foot-points, although 
the stable region is no longer visible. The non-dominant orbits with = 7/11 have moderate distance 
from the Schubart orbit at (f) a = —0.210 and arrive at the boundary of the Schubart region at (g) 
a = —0.230. Since the stable regions are small, they are less effective to the shape of the Schubart region 
than the dominant periodic orbits. The number of arches is five already in (f) a = —0.210. According to 
Paper I, the number of arches becomes five at the totally degenerate case (a '= — 0.208) when the germs 
reg(5, • • •) touch the 0-axis. It is important to note that the value of a such that the dominant periodic 
orbits bifurcate is far from that of total degeneracy. The next dominant orbits with =4/6 exist at 
(h) a = —0.250. In this case, since there is two stable orbits, while their period is three, the number of 
islands is six. This makes, as a result, the separatrix be of hexagram shape. The periodic points with 
^oo = 3/5 still exist at the place marked with itj. 

4 Summary 

We have studied the movement of the periodic points bifurcated from the fixed point (the Schubart 
orbit) and their influence on the structure of the Poincare section for symmetric mass configuration. The 
following is the summary of the present paper. 

• The periodic orbits with v = (n — 2)/n, n being integer, are influential on the structure of the 
Poincare section. There is a rule about the number of orbits for these type orbits. A pair of stable 
and unstable orbits with period n bifurcate for odd n, while two pairs bifurcate for even n. 

• There is a value a — a\/s where the rotation number v at the fixed point is 1/3. Whether a is less 
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(a) ( a ,b) = (-0.150,0) 




(b) (a,b) = (-0.160,0) 






Figure 8: The periodic points and the separatrix on the Poincare section for v 
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(d) (a,b) = (-0.170,0) 




(e) (a,b) = (-0.186,0) 




(f) (a,b) = (-0.210,0) 




Figure 8: (continue) 
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(g) (q,fe) = (-0.230,0) 




(h) (a,b) = (-0.2500,0) 




Figure 8: (continue) 
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or greater than ai/3, as \a — (X1/3I increases, the periodic orbits bifurcate from the Schubart orbit 
one after another. As a result, there are no periodic orbits bifurcated from the Schubart orbit at 
a = a 1/3 . 

• While the periodic points are inside the Schubart region, the germs, which bifurcate from arches 
grow along the separatrix of polygramic shape. The Schubart region takes also polygram shape. 
After the periodic points go out of the Schubart region, the Schubart region becomes polygonal. 
The germs grow along this polygon. 

• The stable periodic points leave the fixed point quickly. They collect germs and sink toward the 
0-axis. This collection of germs results in re-composition of arches. On the other hand, the unstable 
periodic points stay at the vertices of the Schubart region for long time. 

• At the moment of the bifurcation of the dominant periodic orbit for = (n— 2)/n, the number 
of arches is still n — 1. 

Several problems are remained unsolved. First, we could not determine where the stable periodic points 
tend to. One possible answer is to reach the foot-points. Second, the place where the unstable periodic 
points should finally go is also unclear. Third, we have found that two stable orbits bifurcate, when 
(n — 2)/n is divisible by 2. The possibility that three or more (non-dominant) periodic orbits bifurcate is 
unknown. Fourth, we have not considered the periodic orbits which do not bifurcate from the Schubart 
orbit. For example, in Fig. [3(a) a = —0.15, there are four black regions inside arches (one region for one 
arch). We have confirmed that there are two period-2 orbits inside these regions such that one has its 
orbital point in the first and the third arches and the other in the second and the fourth arches. These 
periodic orbits has a different symbol sequence, (21121)°°, from that of the Schubart orbit. 
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